Arbitrary bending of optical solitonic beam regulated by boundary excitations in a doped 

resonant medium 



o 

(N 



Anjan KundtQ 

Theory Division, Saha Institute of Nuclear Physics, 1/AF, Bidhannagar, Kolkata 700064, India 

Tapan NaskaiQ 

Centre for Excellence in Basic Sciences, Health Centre, 
University of Mumbai, Vidhyanagari Campus, Mumbai - 400098, India 
(Dated: March 8, 2013) 

Bending of a shape-invariant optical beam is achieved so far along parabolic or circular curves. We 
propose beam bending along any preassigned curve or surface, controlled by the boundary popu- 
lation inversion of atoms in a doped resonant medium. The optical 2D or 3D beams generated in a 
nonlinear Kerr medium and transmitted through a doped medium represent exact soliton solutions 
of integrable systems. Arbitrary accelerating solitonic beam predicted here, should be realizable ex- 
perimentally and applicable to nonlinear events in other areas like plasma or ocean waves. 
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It is our long-standing dream that one day it would be 
possible to produce localized light beams, that could be 
bent into any preassigned curve with self-acceleration 
and shape-invariance throughout their propagation in a 
medium. Such beams would be able to self -bend around 
an obstacle, making it invisible. However, the first step 
in achieving a self-accelerating beam through theoretical 
and experimental studies has been taken quite late GHH], 
inspired by a pioneering work of Berry et al related to 
accelerating free quantum particles [3]. The beam, typ- 
ically an Airy function solution of linear Schrodinger 
equation, could maintain its parabolic form over a fi- 
nite distance. Subsequently, acceleration along an ar- 
bitrary curve was achieved, though at the cost of non- 
preservation of the shape A beautiful bending in 
circular path is found recently for a shape-preserving 
optical beam on a plane, breathing periodically in time, 
as a Bessel function related exact solution to the lin- 
ear Maxwell equation in vacuum |5|]. A parallel de- 
velopment has also been started in this context involv- 
ing nonlinear media, in which an Airy-like beam solu- 
tion is found to the nonlinear Schrodinger (NLS) equa- 
tion in a Kerr medium Negating the prevailing 
scepticism, that symmetric nonlinearities can not pro- 
duce self-accelerating soliton, such a shape-preserving 
localized static beam solution was obtained (by replac- 
ing time t with optical axis z), oriented along a parabolic 
curve in a two dimensional (2D) plane, together with its 
trivial generalization to three dimensions (3D), which 
subsequently have been verified in real experiments 

Nevertheless, in all such exciting achievements, 
one could obtain shape-preserving beams, which can 



bend basically along two types of curves: circular or 
parabolic, linked to the Airy or Bessel function related 
solutions [l]-[8]. This situation however falls way short 
of our expectation of having a shape-invariant optical 
beam, which could bend to any desired curve, regulated 
by some pre-settings in a nonlinear medium. Such a sta- 
ble self-accelerating beam in 2D or 3D should be able 
to turn around an object of arbitrary shape, turning it 
invisible. Our purpose here is to offer a breakthrough 
by proposing a scenario close to the above expectation, 
where a localized 2D or 3D optical beam could be accel- 
erated in an arbitrary way, controlled by the boundary 
value of the population inversion of atoms in a doped 
resonant medium, through which the shape-preserving 
beam passes and hence can bend around an object with- 
out hitting it. The idea of our construction for bending 
of optical beams follows a fundamentally different path 
than those adopted so far, as briefed above, and is in- 
spired by a series of earlier studies in the nonlinear op- 
tical communication Ir9ll- fl9ll , listed below. 

Soliton based optical communication governed by an 
integrable NLS equation 
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iE z + E tt + 2\E\ z E = Q, 



for the electric field E(t, z) was proposed long back I1C 
where the group velocity dispersion of the field could 
be countered by its self phase modulation in a nonlin- 
ear Kerr medium, resulting to a stable shape-preserving 
pulse propagation in ID [9]. Here the subscripts denote 
partial derivatives, with t as time and z as the coordi- 
nate along the optical axis. The NLS equation (]}, being 
an integrable system, possesses many interesting prop- 
erties like higher conserved quantities, exact N-soliton 
solution, integrable hierarchy and a Lax pair The 
(1 + 1) dimensional NLS equation (l), a successful model 
in nonlinear optics can be derived from the more funda- 
mental Maxwell equation for the electric field in a non- 
linear medium 11211 , where the optical pulse can propa- 
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gate in the form of a stable localized 1-soliton solution: 

E = -2iscch2£e 2ie , (2) 

with£ = r](t—Voz), 8 = —kt+id z r and vanishing electric 
field at the initial moment E(t — s- — oo, z) = 0. Note that 
this conventional soliton moves with a constant veloc- 
ity Vq and modulation frequency uio of the enveloping 
wave, preserving its shape due to constant width and 
the wave vector, all linked to spectral parameter value 
= k + ir]. Though in the context of optical communi- 
cation velocity and frequency parameters: Vo, u>o, repre- 
sent actually their inverses, we continue calling them by 
their usual names, for convenience. 

Another proposal which gained popularity for im- 
proved transmission of nonlinear optical pulses is due to 
self-induced transparency (SIT), produced by the coher- 
ent response of the medium to ultra short pulses |T3l[T^| . 
The governing dynamics is described by the SIT equa- 
tions 

ipt = 2(NE-ui oP ), iN t = E*p-p*E (3) 

together with the coupling to the electric field: iE z = 2p. 
Here p(t, z) is the polarization in the resonant medium 
induced by the propagating electric field E(t, z) and 
— 1 < N(t,z) < +1, is the population inversion (PI) 
of the dopant atoms , ojq being the averaged natural fre- 
quency. The SIT equations is found also to be an inte- 
grable system associated with a Lax pair, admitting ex- 
act localized soliton solutions for all the fields E,p and 
N, involved in the equations with initial conditions 

E ->• 0, p -> and N ->• 1, at t -> -oo. (4) 

The optical pulse described by the electric field can be 
given in the same solitonic form (|2), with changed ar- 
guments as £ = T)(t — v m z), 6 = —kt + uj m z, with con- 
stant velocity v m and modulation frequency u> m depen- 
dent on the constant spectral-related parameters k, p. A 
detailed account of the SIT system can be found in ltl511 . 

A subsequent important idea is to combine these two 
models, e.g. NLS {1} and SIT ((3) equations for propa- 
gation of optical pulses in a coupled NLS-SIT system, 
by transmitting the optical soliton created in a Kerr 
medium further through an Erbium doped nonlinear 
medium flil llTll . Remarkably, this coupled NLS-SIT 
model retains its integrability and could be associated 
with a Lax pair by cleverly combining those of its inte- 
grable subsystems jl6l - [l9ll . The fields E,p, N admit ex- 
act soliton solutions with the optical pulse given again 
in the localized form ((2) with arguments: 

i = r}{t-Vz), 6= -kt+Qz, (5) 

with a constant velocity V = vq + v m and modulation 
frequency = ujq + uj m , given by a superposition of 
the entries from its subsystems. It is crucial to note, that 
the constancy of velocity V and frequency ft is linked 
here to the initial conditions (3). We will see below what 



happens, when the last of these conditions is relaxed. 
Our intention is to use the same governing equations 
but adopt them for the present context of optical beam 
bending on a 2D surface, by formally replacing the time 
variable t with coordinate x along the transverse direc- 
tion: 

iE z +E xx + 2\E\ 2 E = 2p, 
i Px = 2(NE-u p),iN x = E*p-p*E. (6) 

Here p(x, z) = vv* is the polarization in the resonant 
medium induced by the electric field E(x, z) and TV = 
|£| 2 - |^| 2 , -1 < N(x, z) < +1, is the PI profile of 
the effective two-level dopant atoms, with normalized 
wave functions v and v for the ground and the excited 
states, respectively. We recall that the coupled NLS-SIT 
system I0 is an integrable system admitting stable lo- 
calized optical beam in the form of a soliton solution as 
presented in (0 with © (with t replaced by x), for the 
the boundary conditions E — > 0, p — > and N — > No = 
1, at £ — > ±oo (correspond to initial condition (0) for 
x — > t ). Notice that, due the constant PI profile A^ = 1 
at the transverse boundaries, i.e. when the boundary 
atoms along the transverse directions are uniformly ex- 
cited, one obtains soliton solution for the optical beam 
moving with a constant space-velocity V and therefore 
without any bending. 

Arbitrary bending of 2D beam- An intriguing observation 
for this coupled NLS-SIT system ((6) is that, instead of 
choosing the boundary PI as a constant: No = 1, we can 
take the PI at the transverse boundaries as an arbitrary 
function : N = N n {z) along the optical axis z, by suit- 
ably preparing the excitation profile of the dopant atoms 
at the boundaries in the transverse direction. This alter- 
ation in preparing of the set up brings however a signifi- 
cant change in the character of the optical solitonic beam 
propagating through the medium. Since the soliton ve- 
locity is linked directly to the boundary PI configuration 
No(z), the velocity itself becomes a function V(z) of z 
and the exact soliton solution of the system (O under- 
goes a space-acceleration, expressed in the analytic form 

E(x, z) = -2iscch2£e 2?;e , £(x, z) = T](x - V(z)), 
6(x,z) = -kx + n(z). (7) 

Here the variable velocity V(z) and modulation fre- 
quency O(z), given respectively by 

V(z) = v z + v rn f(z), Q(z) = u> z + uj rn f(z), 

where f(z) = [ dz'N (z'), (8) 

are determined by the PI profile, at the transverse 
boundaries: 

N(x ±oo, z) = N (z). (9) 

We call V(z) and fl(z) as the velocity and the fre- 
quency for convenience, though strictly speaking their 
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derivatives V (z) = v a + v m N (z) and Q (z) = u + 
u m N (z) should be the proper entries and be called 
the space-velocity and the space-modulation frequency, 
respectively. Consequently, the nontrivial self space- 
acceleration of the beam: \V (z)\ = -^py , related to the 
curvature R(z), should induce a bending of the optical 
beam. The intensity of the solitonic beam: \E(x, z)\ = 
2sech2£, described by 10 is peaked at £ = = x — V(z) 
and therefore would bend along the curve x = V(z) in 
(x, z) plane, where arbitrary function V(z) is fixed by 
the boundary setting of the PI profile No(z) as in 0. 
Note, that unlike the infinite energy Airy beam, the op- 
tical beam described here by the exact solution 10 is 
exponentially localized and corresponds to a finite en- 
ergy solution of nonlinear equations 10. It is important 
to note for physical applications that, since the curve- 
generating function V(z) is determined by the bound- 
ary PI profile Nq(z), one can tune the form of the curve 
along which the beam would bend, by suitably design- 
ing the PI configuration of the dopant atoms at the trans- 
verse boundaries. 

We present below few examples of such a bending of 
the peak intensity of the solitonic optical beam 10, fol- 
lowing different curves (Fig. 1 a-d) for different choices 
of V(z). The curve-generating arbitrary function V(z) is 
linked to the boundary PI profile Nq{z) as 10 or more 
directly as N (z) = ^(V (z) - v ). 

i) V(z) = z 2 , induces bending of the optical beam 
along a parabolic curve x — z 2 (Fig. la). 

ii) V(z) = (1 — az 2 )^ with a constant a corresponds 
to an elliptic curve x 2 + az 2 = 1, along which the optical 
solitonic beam would bend (Fig. lb). The curve natu- 
rally degenerates to a circle x 2 + z 2 = 1, for a = 1. 

iii) V(z) = z^ 1 yielding zx = 1, describes a hyperbolic 
curve for the bending of the solitonic beam (Fig. lc). 

iv) For the choice V(z) = sinz, the optical beam 
would follow a periodic curve x = sinz (Fig. Id). 

This result demonstrates, that in our setup the bend- 
ing of a nonparaxial optical 2D beam can be achieved 
along any desired planer curve x = V(z), by choosing 
suitably the arbitrary function V(z), linked in turn to the 
boundary PI Nq( z ) 10, through relation 10. Therefore, 
the nature of the curve along which the optical beam 
bends, can be regulated by properly preparing the exci- 
tation profile of the dopant atoms along the boundaries 
in the transverse dimension x. Note, that here we obtain 
a self -accelerating shape-invariant stable solitonic beam 
as an exact solution of integrable equations, seemingly 
contradicting the common belief, that accelerating soli- 
ton can not exist in an integrable system, due to strict 
conservation laws. Resolution of this apparent paradox 
is that, here the acceleration of the soliton is linked to 
a nontrivial boundary function, which serves as an ex- 
ternal source of energy and should also be counted in 
defining the conserved quantities. 

Arbitrary bending of 3D beam- In spite of the success in 
bending of an optical beam along an arbitrary curve 




FIG. 1: Examples of the possible bending of a nonparaxial op- 
tical beam along different curves in x, z plane with optical axis 
z > and transverse dimension x. The figure shows the inten- 
sity \E(x, z)\ of the solitonic beam (7} for different choices of 
the curve-generating function V(z). The nature of the curve 
can be controlled by presetting the population inversion of 
dopant atoms along the boundaries in the transverse direction. 
The panels show the bending of the peak intensity of optical 
beams in the first quadrant of the plane (a) along a parabolic 
curve, (b) along an ellipse, (c) along a hyperbola and (d) along 
a periodic curve: x = sinz. 



on a 2D plane reported above, we intend to general- 
ize this result to include another transverse dimension 
y, since our ultimate interest is to achieve arbitrary beam 
bending in a physical 3D space. This would mean, that 
we could create a stable shape-preserving optical beam, 
that would stretch along the optical axis z with an ar- 
bitrary space acceleration and hence could bend along 
an arbitrary path in a 3D space: (x, y, z). However, 
it is in no way an easy task, since for this one has to 
find first a (2 + 1) dimensional integrable generaliza- 
tion of the (1 + 1) NLS equation 0, allowing stable soli- 
ton solutions, which subsequently has to be coupled to 
the SIT equations in an integrable way, in analogy with 
the integrable NLS-SIT system 10. Note, however that 
the conventional extension of the NLS equation to 
(2 + 1 ) dimensions with cubic nonlinearity will not serve 
our purpose, since it does not allow stable soliton so- 
lutions 12011 . We have to look therefore for an alterna- 
tive higher dimensional NLS equation, that could ac- 
commodate shape-preserving soliton solutions. Fortu- 
nately, such an extension of the NLS equation to (2 + 1) 
dimensions has been found recently l2lt[22ll , by replac- 
ing the standard cubic nonlinearity of the NLS equation 
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by a current like nonlinearity J l E: 

iE z + E yt + 2iJ t E = 0, J* = E*E t — E* t E. (10) 

This 2D NLS equation admits a stable 2D dynamical 
soliton solution for the electric field E(t,y,z) again in 
the same sech form I0 which however is defined now in 
a (2 + 1) dimensional space-time: (y, z, t) with its argu- 
ments given as 

£ = r){t — V\y — V2z), 6 = —kt + u!iy + u!2Z, (11) 

where the velocity components («i, v 2 ) as well as mod- 
ulation frequencies (wi,u)2) are constants and depend 
only on the real and imaginary parts of a fixed spec- 
tral parameter Ai = k + irj . The higher dimen- 
sional NLS equation (flQb is associated with a Lax pair 
({/2(A), V^(A)), containing higher powers in spectral pa- 
rameter A, exhibiting thus an important signature of its 
integrability j2ll[22ll . For confirming the importance of 
this 2D NLS equation in the context of nonlinear optics, 
it can be shown that Eq. (fTOb , similar to the well known 
ID NLS equation, is also derivable from the Maxwell 
equation for the electric field in a medium with nonlin- 
ear polarization H3l . In nonlinear optics the soliton like 
solution i fTTt for the electric field E(t, y, z), would repre- 
sent a shape-preserving stable optical beam transmitted 
through a 2D nonlinear Kerr medium, which can move 
only with a constant velocity v = (vi, V2), as seen from 

However, since our aim is to achieve a stable 3D op- 
tical beam with bending, we have to create an acceler- 
ating soliton solution in place of the constant-velocity 
soliton (TTT l and therefore, in analogy with the coupled 
NLS-SIT system described above, we propose to couple 
the (2 + 1) -dimensional NLS equation (TO) with the SIT 
equations 10 (by replacing t — > x, suitable in our con- 
text) to have a 2DNLS-SIT system as 

iE z + E yx + 2iE(E*E x -EE* X ) = 2p, 

ip x = 2{NE-u p),iN x = E*p-p*E, (12) 

where p(x,y,z) is the polarization in the 3D resonant 
medium induced by the electric field E{x,y,z) and 
— 1 < N(x, y, z) < +1 is the PI profile of the dopant 
atoms as in (6), though generalized here for the 3D 
space. In the set of equations JT21 we have replaced time 
variable t with the transverse dimension x for describ- 
ing nonparaxial light beams along optical axis z with 
two transverse dimensions x, y. We find importantly, 
that the coupled set of 2DNLS-SIT equations fi~2l is an 
integrable system associated with a deformation of the 
2DNLS Lax pair as (U 2 {\), V 3 (A) + V-i(A)), with an ad- 
ditional term V-i(X, N,p,p*), containing the SIT fields: 
N,p,p* and an inverse power of the spectral parame- 
ter: A -1 . Recall that, we have presented above a similar 
NLS-SIT system, but in lower dimensions with the de- 
tails given in fT^,[l9ll . For localized soliton solutions we 
are interested in, the electric and the polarization fields 



should vanish at the boundaries x — > ±00 in the trans- 
verse directions: E(±oo, y, z) —> and p(±oo, y, z) — > 0, 
while the PI with nontrivial boundary values: N(x —> 
±00, y, z) = No, where —1 < N (y, z) < +1, could be, 
in general, an arbitrary function of y and z. We will see 
that, No(y, z) arranged at the transverse boundaries of x 
would be responsible for the arbitrary space acceleration 
of the solitonic beam and hence for its bending along an 
arbitrary surface embedded in a 3D space. 

It is remarkable that, the set of higher dimensional in- 
tegrable equations (12) , similar to the NLS-SIT system 
Ell [III], allow a localized stable soliton solution as 

E(x,y,z) = -2isech2£e 2ie , 
£ = rj(x - viy - v 2 (y, z)),6 = -kx + u>iy + w 2 (y, z), 

(13) 

together with the soliton solutions for the other fields 
p(x, y, z) and N(x, y, z). Note that (T3~3|) is a 3D general- 
ization of 10 and a deformation of (11) . where the veloc- 
ity v 2 (y, z) and the modulation frequency u>2{y, z) com- 
ponents along the z-axis have been changed to include 
y, z-dependence due to more general coordinate depen- 
dant choice for the boundary PI No(y,z), as evident 
from their explicit form 

v 2 {y 1 z) = v 2 z + v 3m f(y 1 z) 1 oj 2 (y,z) =uj 2 z + ui3mf(y 1 z) 

where f(y,z)= f dz'N (y,z') (14) 
Jo 

where the constant parameters v 3m , u>3 m are determined 
through a fixed spectral parameter Ai = k + if]. Since 
the intensity of the exponentially localized solution (13l ) 
is peaked over the surface £(x, y, z) = 0, the localized 
optical beam described by this solution would space- 
accelerate as d zz V2(y,z) along both z and y directions 
and consequently, would self-bend over an arbitrarily 
curved surface x = v\y + V2(y, z). The form of this 3D 
surface is determined by an arbitrary function V2(y, z), 
which in turn is linked to No(y,z) as in (Til l. It is 
therefore possible to bend a 3D optical beam in our ar- 
rangement over any preassigned manifold by regulat- 
ing the setting of the PI function No(y, z), prepared at 
the boundary in the transverse direction x — > ±00. We 
demonstrate in Fig. 2 an example, where the curved sur- 
face along which the 3D beam can bend is described by 
x = y 2 +z 2 , representing a conical surface, spreading out 
with circular projections of increasing radius along the 
transverse direction x > as its central axis, with its ver- 
tex at a: = 0. Such an optical beam can turn around a 3D 
object placed inside the hollow of the cone, making it in- 
visible. In this example, the surface-generating function 
v%(y, z) = y 2 + z 2 — v\y — V2Z, is linked to the boundary 
PI through the relation: N (y, z) = -^{d z v 2 {y, z) - v 2 ). 
The other choices of function V2(y,z) and hence for the 
boundary PI No(y, z) would generate other forms of the 
surface in the 3D space, over which the intensity peak of 
an optical beam would bend. 
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FIG. 2: Conical surface over which a 3D optical solitonic beam 
can bend turning invisible any object placed inside the cone. 
Figure shows the beam intensity \E(x,y, z)\ = sech^x, y,z) 
for a particular choice of function V2(y, z) (see (fT3t'>, where z 
is the optical axis and x, y are transverse dimensions. The 
panels show for clarity different cross-sections of the conical 
surface (for physical interpretation consider only the positive 
quadrants), a) with major vertical section at z = 0.0, b) the 
vertex of the cone at x = 0.0, c) another vertical section of the 
surface at z — 1.3, d) a horizontal section showing a circular 
projection at x = 3.0. 



Based on our findings here we conclude therefore 
that, an optical beam created in a nonlinear Kerr 



medium as a stable localized soliton, can acquire an ar- 
bitrary acceleration while transmitted through an Er- 
bium doped medium, in spite of being an integrable 
system, due to a path-dependant boundary condition. 
This causes therefore bending of an optical beam along 
an arbitrary curve on a 2D plane or over a curved sur- 
face in the 3D space, which can be linked to a nontriv- 
ial population inversion of the effective two-level atoms 
of the doped medium at the boundaries of the trans- 
verse direction. The form of the curve or the surface 
for the bending of the beam therefore, can be controlled 
by carefully exciting the atoms of the doped medium by 
pre-pumping through an external laser. This pumping 
of external energy, in fact, is responsible for the acceler- 
ated motion of the optical solitonic beam. This arrange- 
ment should work both in two and three dimensional 
spaces by coupling the corresponding integrable NLS 
equation with the SIT equations, which would yield ex- 
act finite energy localized beam solutions with self accel- 
eration, inducing arbitrary bending. We hope that this 
possibility, based on the idea used already in nonlinear 
optical communication, for bending arbitrarily a 2D or 
3D optical beam beyond parabolic or a circular curve, 
would be feasible for experimental verification, opening 
up new vistas of practical applications. Moreover, due 
to the importance of the NLS equation and related soli- 
ton solutions in nonlinear phenomena in diverse areas 
in physics like plasma, ocean waves, optical commu- 
nications etc. the possibility of generating accelerating 
solitons in a related nonlinear system should be of sig- 
nificant applicable interest. 
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Supplementary material to Arbitrary bending of 
an optical soliton beam. ...by A. Kundu et al 



I. DERIVATION OF INTEGRABLE ID NLS AND 2D 
NLS EQUATIONS FROM THE MAXWELL EQUATION IN 
A MEDIUM 

A. Maxwell equation for the electric field in a nonlinear 
medium 

Since the propagation of light inside an optical 
medium is guided by Maxwell equation in the absence 
of free charges: 

„ „ „ 1 <9 2 E <9 2 P 

VxVxE= — r^-r-Mo^r, (15) 



c 2 dt 2 



dt 2 



we intend to derive both the one-dimensional (ID) as 
well as the two-dimensional (2D) nonlinear Schrodinger 
(NLS) equation (presented in the main text as Eq. (1) 
and (10)) from the Maxwell equation, to confirm their 
relevance in the context of the nonlinear optics. We 
closely follow the procedure for deriving the equations 
like ID NLS and the derivative NLS as detailed in [1] 
and slightly extend it to accommodate the derivation of 
the 2D NLS. 

Using further VxVxEe V(V • E) - V 2 E = -V 2 E, 
due to V • E = 0, and taking Fourier transform 



E(r, uj - u a ) = / E(r, t) exp — u) Q )t) dt (16) 

J — oo 

equation <fT5t reduces to 

d 2 P 

V 2 E(r,w-cj ) + fcoE(r,w-w ) = A*o^, (17) 

where the polarization is given by 

P = P L + P WL , (18) 
with the linear polarization 

e o f°° 

^ J — oo 

(19) 

and nonlinear polarization 



J — OO 



PjVL = (QX 



with R(t) being the nonlinear response function normal- 
ized as R(t)dt = 1 

Taking the Fourier transform also with respect to y we 
can obtain 



E(p, ft -Ct ,u} - loq) 



E^y.w-woje*^- 00 ^ 

i 

(21) 



where p = (x,z). By adopting the separation of vari- 
able method and the slowly varying envelope approxi- 
mation [1], we obtain 

E(p,Cl-fl ,ui-ui ) = F{x)A(z,n-n ,uj~ujo)e lfioz 

(22) 

Following the details given in [1] one obtains further 
from Eq. ((171 and (|22l the set of equations 



2i(3 A 



e H k 2 - P 2 

/3 2 - Pi 
,(i) 



F 
A 







Po- 



d 2 P J 



NL 



dt 2 



(23) 
(24) 



where e(uj) = 1 + x P ~ an d the scalar 

form of polarization is used with slowly varying optical 
field. Using ll20t therefore, after inverse Fourier trans- 
form in t and y, Eq. d24b can be written as 



dA a A n dA 



+ P< 



01 



dA ip 02 d 2 A 
dy 2 dy 2 



ip2o d 2 A 
2 dt 2 
8 2 A 



+ iP 



dydt 



lJ { 1 + ir ^i) ( A ^ y ' z ^ / R(t')\A(t-t',y,z)\ 2 dt' 



where f3 mn = g^on™ are tne dispersion coefficients and 
Tr = tR(t)dt is the first moment of the nonlinear 
response function. Note that, using \A(t — t',y,z)\ 2 = 
\A(t,y, z)\ 2 — t'dt\A(t,y, z)\ 2 for the slowly varying pulse 
envelope we can rewrite Eq. j25l l as 



(25) 



dA 

dz 



+ -A + /3 



10" 



+P 



01 



dA if3 02 d 2 A 
dy 2 dy 2 



dA t ip 20 d 2 A 

2 dt 2 

■R d " A 

ipiv 



dydt 



r. 



A\ 2 A 



i d 
uiq dt 



(\A\ 2 A) 



d 



T R (A-\A\ 2 



(26) 



where the RHS is simplified and a second order term in- 
volving the ratio has been neglected due to its small- 
ness, as done in [1]. 

From the general form of equation 1 1261 , which is de- 
rived here from the Maxwell equation ((TBl l we can derive 
now ID and 2D NLS equations at different limits, as we 
show below. 



B. Derivation of ID NLS equation 

Considering the propagation of the field A(t, z) only 
along the z-direction, we can drop all terms with ^ in 
d26l l due to the independence of the field on the trans- 
verse dimension. Taking further a loss-free system a = 



7 



with pulses of width T > hps, resulting to very small 
parameters: ^4yr- — > and ^ — )• 0, we can neglect last 
two terms in the RHS of Eq. l l26t to reduce it to 

dA dA i(5 20 d 2 A . 2 

97 + /3l0 9F + — 9^ = ^ (27) 

which shifting to a frame moving with the pulse (re- 
tarded frame) and absorbing the coefficients by scaling 
is transformed to the dimensionless form for the 1 + 1- 
dimensional NLS equation (presented in the main text 
asEq. (1)). 

C. Derivation of 2D NLS equation 

Consider now more general case, when the field 
A(t, y, z), propagating along the optical axis z depends 
additionally on the transverse direction y, with a loss- 
free system a = 0, in the complementary case of ultra- 
local optical pulse with To < lps, together with small 
nonlinear coupling: 7—^0, where we can neglect now 
the first term, but not the last two terms in the RHS of 
Eq. ll26t . Assuming further vanishing values of the dis- 



persion coefficients: f3. l0 -t 0, i = 1,2, f3 02 — > 0, and 
going to a retarded frame of reference t —> t — Pi.qz, the 
general equation j26l l can be reduced to a simpler form 

7T + ^ " 7r£ + ^ ( (~ + iT *) A * At 
oz ay at W^o J 

which by imposing an additional condition Tr = 
reduces further to 

l ^-^dym + l i- A { A m- A ^r)= (29) 

By absorbing the coefficients through scaling of the vari- 
ables, we get finally the required dimensionless form of 
the 2D NLS equation ( presented in the main text as Eq. 
(10)) and coupling it to the self induced transparency 
equations results to the NLS-SIT system. 

This completes the derivation of both 1 + 1 and 1+2 di- 
mensional nonlinear Schrodinger (NLS) equation from 
the Maxwell equation. 



